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ABSTRA.CT 


A theory of the excitation of solar type III bursts is presented. 
Electrons initially unstable to the linear bump-in-tail instability are 
shown to rapidly amplify Langmuir waves to energy densities characteristic 
of strong turbulence. The three-dimensional equations which describe the 
strong coupling (wave-wave) interactions are derived. For parameters 
characteristic of the interplanetary medium the equations reduce to 
one-dimension. In that case the oscillating two stream instability (OTSI) 
is the dominant nonlinear instabxlity. OTSI is stablized through the 
production of nonlinear ion density fluctuations that efficiently scatter 
Langmuir waves out of resonance with the electron beam. An analytical 
model of the electron distribution function is also developed which is 
used to estimate the total energy losses suffered by the electron beam as 
it propagates from the solar corona to 1 A.U. and beyond. 



I. INTRODUCTION 


Among the various forms of active solar radio phenomena^ type III bursts 
are unique in that they may extend throughout the solar corona and inter- 
planetary mediim, to heliocentric distances of 1 A.U. and beyond. These 
bursts are known to be excited by streams of moderately energetic electrons 

which are accelerated either in flares or in active storm regions and re- 

« 

leased along magnetic field Lines that penetrate the high corona. The radio 
emission is the observable by-product of the interactions of the exciting 
stream with the background plasma, and the understanding of these inter- 
actions and of their effects on the propagation of the exciter presents the 
most fundamental problem in type III theory. 

The historically accepted hypothesis is that the exciter constitutes a 
non-thermal feature in the plasma distribution function and is unstable 
against the growth of electrostatic plasma oscillations, or "plasma waves". 
The plasma waves are subsequently converted to electromagnetic waves at 
frequencies either near the electron plasma frequency by scattering from 
the polarization clouds of ions, or near by scattering from each other. 
This hypothesis has been questioned by Kellogg (197^) and by Gurnett and 
Frank (197^), who measured plasma-wave energy densities in situ in type III 
exciters at 1 A.U. and concluded that the turbulence level was several orders 
of magnitude lower than indicated by theory. The latter authors suggested 
that perhaps electromagnetic radiation might result from direct conversion 
of the particle energy and thus be more efficient than the weak- turbulence 
mechanisms implicit in the plasma-wave hypothesis. More recent experiments 
and theory seem to have alleviated these difficulties. First, recent 
observations on board of Helios I and II (Gumett and Anderson, 197^^ 1977) 
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have actually detected wave energy levels in agreement with those predicted 
hy Papadopoulos, Goldstein and Smith (1975j hereinafter referred to as Paper l). 
In addition^ recent theoretical' models involving strong turbulence processes and 
plasma wave collapse* (Smith', Goldstein and Papadopoulos, 197^, hereinafter Paper 
II ) indicate the spatia'l localization of the electrostatic turbulence is such 
that present day experiments have marginal -capability for detecting them. This 
aspect of the problem will be discussed in a companion 'paper (Goldstein, 
Papadopoulos and Smith, 197^) ♦ 

Proceeding on the premise of the plasma-wave hypothesis several authors 
have noted difficulties .in understanding the propagation and dynamics of the 
exciter. Nearly all analyses to date have been done in the context of a one- 
dimensional theory, in which the projection of the particle distribution function 

I 

F (v) onto the direction of the background magnetic field (which guides the beam) 
gives a reduced distribution f (u) such as is 'schematically depicted in Fig. 1. 
The beam-plasma instability leads to the growth of plasma waves with phase 
velocities lying in the resonant region where Sf/du > 0, For homogeneous 
systems quasilinear theory predicts that the distribution function relaxes to a 
plateau (Fig. 2). This led to "S tarrock’s dilemma", after P, Sturrock (1964) 
noted that theoretically the relaxation is expected to occur within a few kilo- 
meters from the injection region, so thatj-type III bursts should never be 
observed. 

This observation led to a series of papers investigating plasma processes 
that could prevent quasilinear relaxation. Two lines -of approach were 
followed. The first one (Group l) was based on homogeneous plasma theory, 
neglecting wave propagation effects but including nonlinear effects within 
the limits of weak turbulence theory (Kadomtsev I 965 ). Prominent among 
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this is nonlinear random phase scattering (i.e.^ induced scattering) on 

electrons, ions or ion sound waves. Kaplan and Tsytovich (I968) considered 

induced scattering off thermal ions (i.e., nonlinear Landau damping on ions) 

as a possible mechanism that could prevent plateau formation. However, as 

was shown in detail (Zheleznyakov and Zaitsev 1970; Melrose 1970a, Smith 

and Fung I 97 I, Heyvaerts and de Genouillac 197^, the time scale for this 

process was much longer than that for plateau formation for type III bursts. 

As noted much earlier by Sturrock (I 965 ); a sufficiently fast time scale for 

nonlinear scattering cauld be achieved, if a sufficient level of ion-acoustic 

turbulence was present. This mechanism will then be similar to the presence 

of a high frequency (u) ) anomalous resistivity (Dawson and Oberman I 965 ) with 

V y, (v __ is the anomalous collision frequency and Y, the instability 

eff b eff b 

growth rate). Following this line Melrose (1970^) proposed that a large 


level of ion acoustic waves could be generated due to the return current in- 
duced by the beam. If the bulk speed of the electrons U = v^ exceeds the 


ion sound speed C = (^)^ v . However, a threshold velocity C requires a 
^ s M e ^ s 

It 

density ratio — > 10 ® which is unrealistically high in the corona. This 
e 

requirement becomes prohibitively severe, if one notes that it is erroneous to 

T 

take C as the threshold for a sr" 0(l) plasma (Papadopoulos 1977)* In this 
® i It _ 

ta TT V sn ilhaii — > lO ^ . Th^ Second line of aDoroach 


case the threshold is U ^ v so that — 

e n 


10 "*^. The second line of approach 


(Group II ) neglected nonlinear scattering effects^ but considered dispersive 


propagation effects of the exciter electrons and of the plasma waves. Within 


this approach Zaitsev^ Mityakov^ and Rapoport (1972) assumed quasilitiear 
relaxation but attempted to take into account the finite spatial extent of 


the beam in a model for its dispersive propagation* Their model is not in 
accord with observations of the exciter evolution at 1 A.U. (kin^ Evans^ and 
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Fainbdrg 1973); their theory has been critized on other grounds by 

, 1 

Smith (1975, 197^). 

Harvey and Aubier (1975) and Aubier (197^) considered the effects of 
wave propagation into a medium of gradually decreasing density. Because the 
frequency of a wave is preserved during propagation^ the wavenumber increases 
in time and the phase velocity decreases; the wave thus passes gradually out 
of resonance with the -particles that excited it. This slowing down, however, 
occurs over distances comparable to that on which the velocity distriWtion 
of the beam must change owing to dispersive propagation, so that the waves 
should continue to be in resonance with some beam particles and only the head 

i. y 

of the beam passes out of resonance with the waves. This fact is not correctly 

T 

accounted for in these papers. Harvey (197^a^, b) constructed a model for 
the exciter evolution In both time and space, including propagation effects. 
Although this model is detailed and elegant, it encounters difficulties 
which may limit its self-consistency by neglecting to account for processes 
such as pitch-angle scattering and resonance broadening, which should be 
incorporated, at least phenomenologically, in such a model. 

In Paper I it was noted that the plasma wave levels expected from all the 
above theories were such that nonlinear, coherent effects associated with 
strong turbulence theory, such as parametric instabilities, plasma wave collapse 

r 

f 

and density cavities will dominate, thereby invalidating their calculations. 

The main nonlinear effects which can be experienced by the electron plasma waves 
were included in the Appendix of Paper I. These equations reduce to a non- 
linear Shroedinger equation when ion inertia is neglected, which was first studied 
by Zakharov (1972). To simplify the formalism it was noted that for type III burst 
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parameters the nonlinear equations of Paper I, would be equivalent to the 
parametric oscillating two stream instability. Notice that as shown in 
Manheimer and Papadopoulos (1975) soliton formation and wave collapse are 
nothing more than the description in configuration space of well-known 
parametric instabilities discussed in Paper I. These instabilities (or 
equivalently the spatial collapse of the waves) act to transfer waves in 
fc-space, from the linearly unstable beam- resonant region to lower, non- 
resonant phase velocities. The OTSI produces S3m3metric spectra of plasma 
waves at both positive and negative phase velocities in the non-resonant 
region, as well as purely growing (i.e,, aperiodic) ion density waves. The 
object of Paper I was to demonstrate that the OTSI may play an important 
role in the interaction of the type III exciter with the ambient plasma. 

With that aim. Paper I considered the temporal instability of a uniform 
beam, with the beam-resonant waves treated as constituting a monochromatic 
pump with = 0 . The use of = 0 (dipole approximation) allowed us to 
concentrate on spectral transfers that, contrary to the weak turbulence 
theory, go towards shorter wavelengths similar to wave collapse. 

In this paper we extend the analysis of Paper I, In Sec. II we briefly 
describe the quasilinear theory of the injection of an electron beam into 
a half-infinite space. The levels of Langmuir turbulence produced in that 
situation are found to be characteristic of strong turbulence; and thus use 
of the quasilinear theory, which ignores all strong-coupling (wave-wave 
interactions), is generally invalid when the dynamical evolution of the beam 
at a particular spatial location is considered. Thus, an essential differ- 
ence between this work and those of Group I is that we consider only strong 
turbulence processes which are ignored in previous analyses. Because the 

ORIGINAL PAG® IS , , 
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physical situation is dominated by strong turbulence, we do not consider 

scattering by weak 'turbulence . Moreover, in contrast to the papers o£ Group 

we do not consider wave propagation effects; onr theory is strictly local* 

The justification for this procedure is that the time scale At of the spectral 

2 

evolution is such that v* At « L , where—v -v /-v, is the group velocity of 

.. g . ti-’ g e h 

the plasma waves and L is the scale length of the solar wind density gradient 
r * n 

(Brejzman and Ryutov 197^). 

In addition our computations do not consider effects of possible density 

inhomogeneities created by the beam-plasma interaction itself, except in-so- 

* 

far-as the strong turbulence modifies the linear Bohm-Gross dispersion 
relation. 

As in our previous work, the theory is one-dimensional. The latter 

feature has been critized recently by Bardwell and Goldman (1976), who argue 

that the interaction is significantly modified by consideration of modulational 

instabilities in three dimensions. We show in Sec. Ill that when the three 

dimensional effects are properly analyzed the spectrum of pump waves excited 

by the electron beam is indeed predominately one dimensional as we previously 

assumed. Furthermore/' the'daughter waves resulting from the OTSI are shown 

to propagate predominately parallel and anti-parallel to the direction of the 

beam. Finally, the decay instability is shown to have a higher threshold 

than the OTSI, at least for parameters characteristic of the interplanetary 

medium. Sec. Ill contains also a derivation of the corrections to the Bohm- 

Gross dispersion relation of linear Langmuir waves which are of great xm- ; 

portance in the presence of nonlinear ion density fluctuations, such as are 

excited by the OTSI or collapse. This dispersion relation is shown to 

eventually lower the initial threshold for the OTSI; which in turn enhances 
' . . 1 i ' 
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the scattering of the pump waves to low phase velocities out of resonance 
with the electron beam. The final stabilization of the OTSI is accomplished 
by an anomalous resistivity that is also excited by the ion fluctuations in 
a fashion similar to Sturrock's. In the following section, we develop an 
analytical model that describes the total energy loss of an electron beam 
as it propagates from the corona to 1 A.U. 

In the accompanying paper (Goldstein, et al . 1978) the theoretical * 
equations derived here are solved numerically, using parameters appropriate 
to interplanetary conditions between 0.1 and 1 A.U. Conq>arisons between the 
predictions of our theory and observations as well as the radiation calcula- 
tions appear in this paper that constitutes the second part of our work. 


ORIGINAL PAGE IS 
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II. LIMITATIONS OF QUASILINEAR THEORY 


The fact that the type III exciter is not a uniform infinite beam, 
but rather evolves continuously in both coordinate and velocity space, has 
important consequences for the nature of the beam-plasma interaction. Re- 
gardless of the details of the beam evolution, two consequences in partic- 
ular may be deduced. The first is that at any point in space, the plasma- 
wave energy density may exceed by a large factor the limit imposed by quasi- 
linear theory for the uniform case. Notwithstanding this fact, however, 
the head of the beam does not suffer significant quasilinear relaxation. 

To better understand these assertions, consider first the simple model 

Jt 

of a beam such as in Fig. 1 injected into a half-space (Fig. Ja). In 

this situation the phenomenon of "oscillation pile-up" occurs (Tsytovich 

1970 j pp. 177)* The unstable waves are convected away from the injection 

point at the group velocity v , while the beam particles move at v. » v . 

g b g 

Conversely, the continuous injection 'of "fresh" particles into the plasma 
may build up the local wave energy density beyond the limit of the uniform 
case, in which it obviously cannot exceed the energy density of the beam. 
This buildup is asymptotically limited by the competition between convection 
and instability. The asymptotic energy density is given (Tsytovich 1970 ) 


W 


( 2 . 1 ) 
v2 


= I? “ “b'nH 

2 

which represents an enhancement of the uniform limit by the factor (v. /v ) . 

D 6 

Physically the reason for this enhancement is the fact that in a steady 

state injection the important quantities to be balanced are the energy fluxes 
1 3 

Ci.e., v ~ level is reached within a distance 



from the injection point, where 


(\ 

% V^’'i 


0) 


(2.5) 


is the coarse-grained average of the heam-piasina instability growth rate.! 


The distance X is limited at any point x to 


X (x) < 


/ Av, \2 /v \ h 


(2.4) 


where is the electron Debye length. For ~ c/3^ ^^b^^b ~ gives 

values for X of less than about 60 km at the 100 MHz level. Even at 1 A.U. 

X < 10~^ A.U. Because these values are much smaller than the corresponding 
length scales of the ambient density gradient, the injection model of oscilla- 
tion pileup is at least qualitatively appropriate. Then Eq. (2.1) implies 


that values 


W 


H ■ y .— Pi tl f— 

n T \T / 

e e \ e / 


(2.5) 


may be attained, where Tl s n^/n^,E^ is the nominal beam energy. At 1 A.U., 
the observations of Lin, Evans, and Fainberg (1975) sod of Frank and Gumett 
( 1972 ) imply T] ^ 10~®, while eq. (2.5)j will give W » 1. The significance 
of this fact is that such high energy densities in the resonant wave spectrum 
are sufficient to drive fixed-phase instabilities even with broadband pump 
spectra. 


Notice that the presence of plasma waves with energy W, will produce 

6n W 

in an isothermal plasma a density depression — = due to pressure balance 
(nT + W = const.) Since = ^ (l + §kX ), the nonlinear plasma frequency 
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will be \ = + l - I ~ 




1 ^ ^ 

2 nT /. 


Therefore-, 


'when W « , strong turbulence phenomena such as OTSI or collapse 

will become important. Since for type III bursts (kX^)^ 10 ^ at 

1 A.U., this level which constitutes the threshold for our processes is 


easily achieved. A more detailed derivation of the nonlinear modifications to 


the Bohm-Gross dispersion relation is deferred to Sec. HI. 
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III. INSTABILITY THEORY 


In Papers I and 11^ we utilized a simple derivation of the OTSI growth rate 
which assumed a monochromatic, one dimensional pump wave spectrum. We give here a 
detailed derivation of the general dispersion relation for modulation and collapse 
instabilities driven by an arbitrary spectrum of plasma waves, For the para- 
meter range of interest for solar type III bursts, we show that one is justi- 
fied in approximating a general three dimensional pump spectrum with a one- 
dimensional one. Furthermore, the fastest growing daughter Langmuir waves ^ 
are shown to also be one dimensional; aligned approximately parallel and 
antiparallel to the mean magnetic field. The general OTSI dispersion re- 
lation thus reduces to the relatively simple form used in Paper II; and so 
we again find that the OTSI provides a mechanism for stabilizing the bursts 
against catastrophic energy losses that would otherwise prevent them from 
travelling out to 1 A.U. 

We conclude this theoretical development with a discussion of additional 
nonlinear effects that are important to the eventual stabilization of the - 
OTSI, The first is that the Bohm-Gross dispersion relation is changed by 
the density fluctuations amplified through the OTSI or their subsequent 
decay to ion acoustic waves . Then these density fluctuations produce an 
anomalous resistivity that enhances the tendency of plasma waves to scatter 
to shorter wave lengths. This in turn hastens the decoupling of the electron 
beam from the Langmuir radiation. 

Note that the effects described here are completely nonlinear and our linear- 
ization corresponds to determining the time scale of the events. 

% 
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These diverse phenomena can be described in terms of a two-fluid hydro- 
dynamic model of the plasma in which the motions of the plasma are separated* 
into fast plasma oscillations (a) (O ) and quasineutrai slow sound oscilla- 

* r C 

tions. We take the ion density as 

= Hq + 6n^ ^ (5*^^ 

and the electron density as 

n = n, + 6n 6n /n. « 1 (3.2) 

e 1 e e r / 

For the high frequency fields the linearized electron equations are 


with 


S6n 
e 

9t 


- - V» (6nj6v^) 


1 “e 


a Y 

n ^ 6v + ^ V6n = -~n e-vn 6v 

oot-e m e - 


e 

m o 


e o -e 


( 3 . 5 ) 

(3.^^) 


Ve = - 4TTe6n^ (5.5) 

and where is bhe polytropic index and is a phenomenological damping 
decrement. With the right hand sides equal to zero, these equations reduce 
to the usual second order equations that describe Langmuir waves. 

The equations for the ions are similar except that the ions experience 
a low frequency pondermotive force, F, caused by the gradient of the pressure 
of the Langmuir waves. If we define P = to be the total plasma 

pressure, then 

F = . ® (£ ) ^ <gg) 

2m ™ - — 

e 

where < > denotes an average over the fast time scale. The linearized ion 
momentum equation can now be written as 
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dv 

nM^=-VP + F - v.itiMv. (5,7) 

dt - 1 -i I / 

2 2 

Using an equation of state of the form P = Me n with c = (y T -f Y.T. )/M 

s s e e XI. 

in conjunction, with the continuity equation^ the ion fluctuations are found 
to satisfy 


+V A 


Under the restriction that (kX 


(5*^) can be combined, yielding 


2 2 

C"" 6n. 

s 1 

« 1 and 



V<ee) 


(5.6) 


^ « 1, equations (5-5) and 
e 


Q 2 e at m e 


U) 

e ^ dn. e 

n 1 - 
o 


(3.9) 


Equations (5*6) and (5*9) ^xe the fundamental equations that describe the 
coupling of high frequency plasma waves to low frequency density fluctuations. 
The equations can be simplified by introducing the slowly varying quantity 
E(x,t) defined by 


e (x,t) =2 




^t)e e -I- E (x^t)e e 


i03 t 


Iti terms of E(x^t) equations (5 *8) (5*9) hecome 

U) 


a 

e 


E = 6n.E 
— 2n X— 


and 






( 5 . 10 ) 


(3.11) 


(5.12) 


Equations ( 5 .II) and (5-12) are the basic equations that describe either the 
OTSI or the formation of soli tons and cavitons (Manheimer and Papadopouios 1973). 
In this section we are primarily interested in, deriving the growth rate of the 
OTSI and comparing it with similar growth rates for other modulation 
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instabilities. The spatial structures are discussed in Goldstein et al. (I 978 ) 


Upon defining the Fourier transform as 


{8n(k,t), E(k,t)} = — ~ Jdx e"^-*? | 6n(x,t), E(x,t)}, ( 3 - 15 ) 


(2 IT) 

equations (3.II) and (5.12) become 


(5.14 


and 


S .. .2 P' 

— o + ^ . — + k C 
5t^ ^ at 


. 2 * 

5 (k,t) = /dk'E(kO*E (k-k') (5.13) 

i6TTMy “ “ 


A coiqjlete solution of these equations including all the nonlinearities 
is beyond the scope of this discussion, for here we are most interested in 
studying the linear stability properties of a Spectrum of Langmuir waves that 
is excited by an electron beam in order to determine which of the various 
modulational instabilities is likely to dominate under plasma conditions 
characteristic of type III bursts. We assxime that a pump spectrum of Langmuir 
waves exists between < k < k^. These waves are taken to be zero-order 
quantities, while waves outside that region are the small perturbations 
6E(k,t) and fin(k,t). From equations (3.14) and (3.15) we find 


" '"ek 


- I t«g(kX^)^j SE(k,t) = ^ Jdk'6n(k-k',t)E(k' 


,t) ( 5 . 16 ) 


and 

r 


a" 




^ + k^C^ 


-k ^ r 

6^(k,t) = E(k;t) • 6E (k-k,t) + C.C. 


(5. IT) 


We now define the propagators L^(k,t) and Lj^(k, t) through the relations 
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( 5 - l 8 a ) 


U) P 

Lg(k,t) 6E(k^t) = ^ 6a{k-k‘',t)E(k',t). 


I^(k,t) 6M(k,t) = Jdk [E(fc%t).5E*(k-k",t) + C.G.] (5.l8b) 


16 ttm 


From ( 5 . 16 ) - (3.18) we obtain the dispersion relation 

2 

[”k •(k-k^)l ^ 5j5 


k^ou 


I^(k,uj) = dkME(k')l‘ 


52 lTMn V 
o 


Tk ' ( k - k ') ^ 

ifqW ” 


(3.19) 


-1 


where is the inverse of the propagator L^, the integration is over the 

volume in k-space initially occupied by the unstable waves, and U)(k) = 

2 

5/2 «>g(k^-^) . From the definitions of and (3*19) can be rewritten as 

2 


r 0^1 r - 

L-a>(a>fiv . ) + k c J = - JdklE(k01 X 


X 


[k'-(k-k')]^ P o n-1 

[- i * i (b-S')\“e ] ^ 




-1 


(5.20) 

The theory as originally presented in Papers I and II was one dimensional. 
Bardwell and Goldman (I976) have maintained that a three dimensional stability 
analysis is essential. In asking whether a one dimensional treatment is 
allowed, there are actually three separate questions that should be addressed: 
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l) Is the spec.trum of waves excited hy-the electron beam predominately one 
dimensional? 2) Are the daughter waves produced by the parametric instabil- 
ities predominately one dimensional? 3) What is the role of the decay 
instability? 

We discuss each question separately below, 

l) The growth rate of the beam plasma instability in the presence of 
a magnetic field and for arbitrary direction is given by li 


n^ tu-k*V^ 

Vn. = /tt exp(-X)l^U) ■ X exp 


n 


/ (JU-k-V, >,2 

(^) 


( 3 . 21 ) 


k AV 

where X p** and we assume that the beam velocity spread is isotropic 


e 

given by Eq. (5.22) allows us to calculate the angular spread of the 

spectrum, since Y^(k^)+0 only if X « 1. Therefore, the angle 9 will be 
given by 


sin0 




U) 

e 


1 ^ 

1 _e 

3 u) 

A 


( 5 . 22 ) 


n 

Since — ^ 10 , the angular spread of the instability spectrum is leas 

e / n ' \ 

than 1°. Namely the magnetic field despite its low value ^ 10 forces 

e 

the system to behave one dimensionally. 

2) From the general dispersion relation, equation (5«20), we can now 
assume that the pump spectrum W(k'*) is one dimensional to a high degree of 
approximation. Therefore, in spherical coordinates we write the spectrum of 
pump^ waves as v - . ^ . 


W(k') 


W Ak^6(cos9-l) Ak 

« 1 

2tt\'^ (k -k )^+(%) 


(5.25) 
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where is the width of the pump wave spectrum centered at and now is 
assumed to be aligned along the beam direction. For convenience we will use 
a dimensionless notation in which (kl (u (u/w , v-+v/o) , € = m/M, 



T. 

and T = 1 + I ^ . 


i 

From equations (3*^) and (3*25 ), we find on integrating over the 
resonance that 


)-ek^3 + 


r ek\F(k,k 4) 


^ k^-C (uH-iv^ )+i5 1 kl Ak^cosiir-3kk^cosi}i]^ 


= 0 


(3.2h) 


where i|f is the angle between k and the magnetic field direction, F(i]j) 

“ 2 2 

— cost'll for k » k and F(ilt) — 1 - k sin t k « k . In the 

o ^2 o 

o 

limit (Ak^/k^) 0, equation (5.24) reduces to the dispersion relation found 

by Papadopoulos (1975) and Bardwell and Goldman (1976). It is interesting to 
note that the finite width of the spectrum enters like an electron dissipation 
C i .e ♦ lii + 03 4* i('^g + 3kAk^cos'|f )1 . Physically, this is due to a decorrela- 

tion of the waves produced hy the finite bandwidth. For parameters of interest 


xn our case 


, W < S and the ion response can be neglected Ctu(o>fiv.) < €k^3 , 


so that the dispersion relation becomes 


4 2 ^ ^ 

k - [oH*5ikAk cos'jf-3kk cos'!;] - ^ ° 

O 0 ' i 


(5.23) 


Therefore, 


0 ) 


= 5kk cos'll - 5 IklAk cos'it+'5 i [W - 5k (3.26) 

O O ^ O -L 
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The maximum growth rate as a function of angle is given by S(lirfD)/5’li = Oj or 


- 31 k] sin'Jt 


/- r 


WF(t) u 
^ 5k^ 


’ '•-s 




(5.27) 


which can be satisfied only if ijr = 0 for both k < k^. Neglecting ion inertia 

r p 

implies that k^ < m/M as w^lj. as < ra/K, both easily satisfied for type 
bursts near 1 AiU, We therefore conclude that the one dimensional treatment 
of parametric instabilities is easily justified for the case under considers-^ 
tion. 


5) The final question which we want to investigate is whether the decay 

-5 

instability is important when — 10 . This is a three wave interaction 

'whose threshold is found by setting F (\[t) = 1 and cos i]; = 1 in equation (5*2h). 
The result is 


^Lnr ^ g !l_ 
decay ekT 


(3lk|Ak X 
o 


For T It. ^ I, v./ekT ^ 2 and 

e 3. i decay 


50 |klAk » 50 k Ak = 50 k‘^ Ak /k 

o oo o oo 


^ 50 (l/ 5 )k^ lOk^ » 10~^. 

Thus we have not exceeded threshold for the decay instability. We con- 

I 

elude then^ that for the case under consideration our assumptions are very 


well justified. This allows us to proceed as in Papers I and II without 


the necessary encumbrance of using the three dimensional dispersion relation. 


If w'e define the frequency shift between the pump waye with wave number 

k and the daughter wave with wave number k^ as 6 = - ^ek^ where both 

o 

tu and / satisfy the Bohm-Gross dispersion relation, and 6 <0 for the 

CiX 

o 

OTSI; then the growth rate for the OTSI can he found from (5*25) to be 


approximately given by 
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OTS ’ 0 


(k ,k’') = -l/2 (u)f + 6^) + 


26 ( 1 ) . ^ 

p efc 


X 

s 


, 2 in 

where (o = — 
K 


(k X ) 2 

U)^ 

^2 e 


(5.28) 


l+(k%)' 

From (5*28) the threshold for the OTSI, in the absence of collisions is given 


by 


= - 26(1) , //u) 

T ek e 


(3.29) 


The amplification of aperiodic sound waves by the OTSI eventually changes 
the real part of the plasma dispersion relation of the Langmuir waves, and 
this will in turn modify the threshold condition (3.29). This effect has 
been studied in the dipole approximation by Kaw, Lin, and Dawson (1973) 
albeit in a different context. We present here a generalization of their 
analysis appropriate for a finite wave number pump. Starting with equation 
(3.9) 01^® assumes the existence of an ion density fluctuation given by 


6n. = n (1 + 5 cos(k.x)) 

X 0 1 


(3 .50) 


which, with = 0, gives 




3^6 


— — ^ + u) (l+^cosk.x) S — ' ' ' — n — 0 

3? ® ^ ® 3x2 


(5.51) 


Again we have restricted our attention to one dimension and have used an 
adiabatic poly tropic index. Taking the Fourier space and time transforms of 


(5.31) gives' . 

t t ; » 


2 2 

V» 

e ^ m 


2„ 


^jj(k) + 5 [6p(k+k^)+€p(k-k^^] = 0 


(3 ..52) 


Because c^ouples to and €^(k^k^) then couples to €^(2k^k^)^ 

etc*; equation (5»32) in fact represents an infinite system of equations. 
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For simplicity^ we follow Kaw^ Lin, and Dawson (1973) restrict ourselves to 


coupling only between k and dispersion relation can then be 


written as 




e 


10 , 2/2 




e 






( 3 . 35 ) 


Note that for k=0 we recover the results of Kaw^ Lin and Dawson which 
they derived by solving the (Mathieu) equation for As it is not our 

purpose in this paper to present a detailed discussion of the scattering of 
Langmuir waves from ion fluctuations^ we give only an approximate solution 
of the cubic equation for ( 3 . 53 )j and proceed by assuming that a wave at 
wave number k couples only to ( l.e. , to a wave of lower phase velocity) 

rather than to kfk^. This approximation reduces (5*33) to a quadratic equa- 

p 

tioh in and' introduces an error of order two in the nonlinear frequency 
shift which we are about to derive - as can be confirmed by returning to the 
dipole approximation, k= 0 , and then solving ( 3 * 33 ) both exactly and with the 
approximation introduced above. 

Equation (5.53) then reduces to 

_e 

2 ~ • ■''-- q ' ■ ■'<- '"i’e 


(ju“ 2 

® ^ 1 + 3(k\ ) + 5/2 (k.X^) (k +k) \ 

cu- ® i 6 X € 
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. 2 2 P ’2 j 1/2 

^5/2 T/ 




From ( 5 , 29 ) it is clear that if one chooses the "plus" sign in (3.5^) to go 
with the pump wave and the "minus" sign for the daughter wave, any significant 
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excitation of nonthermal ion fluctuations by the OTSI will tend to decrease 
l^l and thus lower the threshold for the OTSI. 

A further effect which we consider is that of anomalous resistivity 
(v. Dawson I968). In principle, this should be contained in the analysis 
described above. However, to exhibit that, it would be necessary to further 
generalize the dispersion relaxation (3 •53) to include coupling to many 
modes 1), Because that is beyond the scope of the present discussion, 

we have chosen to utilize equation (3 •5^) to illustrate only the physical 
phenomenon of a reduced threshold for the OTSI in the presence of a nonthermal 
level of density fluctuations. We describe the subsequent dissipation of the 
long-wave length waves to shorter-wavelengths in terms of the high-frequency 
anomalous resistivity (Dawson and Oberraan I963), which is excited near m 

e 

by the ion fluctuations. In this process, the electron energy originally 
carried by long wavelength Langmuir waves is conserved (because the ions are 
treated as fixed), but cascades to shorter wavelengths as the electrons 
scatter from the correlated ions. In Fourier space the spectrum of Langmuir 
waves will contain spectral components at high wavenumber. [Any new com- 
ponents appearing at small wavenumber would ultimately be rescattered, either 
by the OTSI or by further anomalous dissipation back to higher wavenumbers, 
and so we will not consider this possibility in our subsequent analysis.] 

A complete discussion of plasma anomalous impedance excited near 

by oscillating electric fields is contained in Dawson and Oberman ( 19 ^ 3 ) 

< 

and Dawson (I968); here we only outline the derivation. 

In a plasma of isotropically distributed ions the effective electron- 
ion collision time T^, is related to the scalar in^edance by 

Z(m) = (l-i/m ) (3.55) 

(JJ 

e 
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[S(k.)-S^(k.)] 


(J.38) 


where S^(k) is the thermal noise level of the ion spectrum. 

A complete description of this non- stationary turbulence should include 
a mechanism for the damping of the aperiodic ion fluctuation in the presence 
of the Langmuir solitons when the electrostatic fields have fallen below 
threshold for the OTSI, A theory of this highly nonlinear situation has yet 
to be developed, though various aspects of the problem have been studied. 

For exa3i5>le, Yu and Spatschek (1976) have considered the damping of ion 
' acoust^ic waves in a medium containing strong Langmuir turbulence. In their 
analysis, the ion waves are assumed to satisfy a normal mode dispersion 
relation; whereas, the ion -fluctuations excited by the OTSI are not normal 
modes, but must have frequencies «•' < In fact, we cannot make use of 

the Yu and Spatschek result in the problem at hand. 

- — However j-a-detaiied knowledge of the ion fluctuation damping mechanism 
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is not really ^necessary in the present discussion so we have adopted the 
following approximations j which we feel do no violence to the basic physics. 


For T > T . we find 
e X 



If T « T. 
e X 



e 


The electron oscillations are also damped by the thermal plasma, 
the solar wind the electron distribution function is characterized by 
thermal component at a 1.15 x 10^ K and a suprathermal tail (0 ^ 
which we have approximated by (Linj Anderson and Cline 19T2). 


1.07 X 10“^^ 11^ 

f = 

1.02 X 10*”^ 

-it 

where « 5*5 x lO” /n and llg = 2,5 x 10”°/n. 

1 

Defining 


1.01 ^ IpI ^ 0.088 

IpI s 0.o88 


f (P ) + £(P ) 

w 

f 

w ^ 


P > 0 
P < 0 


(3.39) 


In 

a 

0 . 01 ) 


Oih-l) 


(3.^2) 
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We can include Landau dancing by aolar wind electrons along with linear 
growth of Langmuir waves excited by the electron beam and reabsorption caused 
by evolution of the beam to lower velocities, by evaluating 



P = /kc 


(3A5) 


The transfer rates Yt^ yrvr<if and v enable us to describe the 

It NL S 

evolution of the turbulence sprectnm by a set of rate equations similar to 
those used in Paper II. As in that discussion we denote by P(k) and A(k) 
the plastnon and phonon number densities which correspond to the energy 
densities W(k) and S(k): 


W(k) = P(k)iio) , /n T 
ek e e 


S(k) = A(k)Hu) /n T 


( 3 . 44 ) 


The rate equations are then 
dA 




dt 


I " Vg(h)|A(k) 


(5.45) 


^ = YL(k)P(k) - S Yo^s(k,k')A(k') - P(k) S YjjL(k') + 


fc >k 




k<k 


( 3 . 46 ) ^ 


df * Yj^(-k)P(-k) - S YQjg(-k^k )A(k ) - P(-k) S Y«t (“^^ ) + 


k >k 


k >k 


NL' 




k^k 


(3.47) 
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In writing (3 .45) - (3.4-7) we have normalized all the transfer rates to 
and have written A(k) = A(+k) + A(-k) so that all k and k are assumed 
positive. A detailed discussion of the results of the integration and their 
relationship to type III hursts observed near 1 A.U. is given in the 
companion paper (Goldstein, et al . l^jQ) . 

However, before concluding this description of the basic ideas which 
comprise our theory of type III bursts,, it is necessary to discuss the 

} 

analytical model which we have used to describe the temporal evolution of 
the electron exciter distribution function. Once this is accon 5 >lished, it 
is then possible to estimate the total amount of energy lost by the electron 
beam as it propagates to 1 A.U. from the sun. These are the subjects of 
the next two sections. 
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IV. evolution of the exciter distribution function 


The beam distribution used iii our calculations has been modelled semi- 
empirically from in situ particle observations in type III sources at 1 A.U. 
(Lin 1973j Lin, Evans, and Fainberg 1973)* Such a model cannot be developed 
from a priori argxments because the propagation of the exciter in the corona 
and interplanetary medium is affected by magnetic field pitch-angle scatter- 
ing, wave-particle interactions, and the inherent dispersion of the injected 
distribution. The effects of these phenomena are manifested on disparate 
length and time scales. Thus, any attempt to model the beam evolution must 
be somewhat heuristic. 

Lin observed the differential energy flux of the beam electrons rather 
than the velocity distribution. His detector accepted only 

large pitch-angle particles (roughly 60 ° to 120°) because the axis of the 
acceptance cone was oriented perpendicular to the ecliptic plane, while 
the magnetic field was predominatly in the ecliptic plane. A second detector, 
with less energy resolution, provided pitch-angle distributions with resolu- 
tion of about 50°. Time sequences of both energy flux distributions and 
pitch-angle distributions during several events were provided to us by Lin 
(private communications). While these distributions are difficult to model 
analytically owing to lack of resolution, they do indicate, however, that 
marked anisotropy occurs early in the event and persists for some time. 

Thus, we have interpreted the differential energy flux distributions as 


one-dimensional velocity distributions. Within this context, the parameters 
to be deduced for the beam are and Av^ all as functions i^f time. Me 

sketch the method briefly here, a more detailed description, including a dis- 
cussion of the interpretation of the spectra as one -dimensional distributions, 
has been given in Taper ~H-. ■ ~ - 
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Lxn^ Evans ^ and Faitiberg (1975) showed that the first particles of 

I 

any given velocity to reach their detector all appeared to have traversed 
the same distance L after injection at a connaon time t^. Thus at time t, 
the slowest particles have speeds given by 






(4.1) 


The peak of the distribution is at v.^^ = P^c + Av.j^. Empirically, we find ‘that 
v^ can be approximated for a wide range of events by 



= PoC ®(P^) 

(4.2) 

and so 




Av, , 



b _ 1 1 

’'b w 

(4.5) 

where ®(P ) 
o 

= 1.4 (140.075?^ + 0.085S^^) (0.3/P 

(4.4) 


To model T1(t), we use the fact that the fully-evolved beam is a power 
law distribution in energy, f(E) ~ E ^ . The corresponding one-dimensional 

distribution is 

f^(p) - A(P)p-^ (1..5) 

where the normalization constant A(p) is defined through 

00 

J dpf^(p) = n^(P ^ P), (^.6) 

? 

P is defined in5>licitly through equation 
and is determined for various events by integration of the observed 
spectra over energy. This latter procedure involves assuming a form for 
the pitch-angle distribution, but a mlnlmian density is attained by integrating 



the spectra observed, by bin with the assumption that it is isotropic. The 
effect of this is to tend to underestimate the growth rate of the linear 
instability^ we shall consider the consequences of this below. 

We asstirae that the peak of the distribution, at P = P^(t), is at the 
fully-evolved value implied by equations’ '‘5") (’4,6), The distribution - 

for ^ P ^ modelled by a tamp, while for p ^ P^ it is given 

by equation (4,5)* Thus, 




(4.7) 


(P-P ) 

fl(p) = f (P.) ^ ^ ^ 

i ° ^ (P.-P ) ° ^ 

P O 

The evolution model of equations (4,l) - (4.4) was developed from obser- 
vations at 1 A..U. in order that the ntomerlcal results of our calculations 
may be compared to the corresponding radio observations. In this context, 
we may then ask to what extent the model permits such comparisons and \diat 
are the uncertainties involved? There are essentially three points to con- 
sider: 


(i) How valid is the direct interpretation of the differential 
energy fltix dj/dE as a one-dimensional distribution 
function f(u) mdj/dE? 

(ii‘) For calculations using the measured flux distribution para- 
meters of a particular event, how accurate is the determination 
of the density normalization in equation (4.6)? 

(iii) IJhat is the validity and what are the consequences of assuming 
a one-dimensional ran^ for (p), as in equation (4.7)? 

Wi th re gard to point (i), our assumption is equivalent to assuming that 
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the distribution is dominated by small pitch-angle particles at any given 

energy. Assuming that the differential energy flux distribution is separ- 

2 

able into the form (for E > E = inv /2) 

b o 


= AE- 

dEdn 


G(y.)H(E-E^) 


( 1 ^. 8 ) 


where “"p. = cos a, gyrotropy is assumed, and H(x) is the Heaviside step function; 
f(u) is given by ^ 

£(u) °= I dw ^ g(u/v) 


where 


X = 


u (u ^ V, ) 
D 

''b ''b’ 


Talcing f(u) mdj/dE implies that <?(P<) is sharply peaked near p = 1. If G(p) 
is expanded in Legendre polynominals, one can show that the peak of f(u) is 
preserved at u = v^, while the slope for u < v^ depends on the relative 
weights of the various terms in the expansion. For example, if G(p-) = 

(1 + p)/2, then f(u) is as shown in Figure (4a), while if G(P') is dominated 
by larger-order (n > l) odd pol 3 niominals, f(u) may be as depicted schematically 
in Figure 4b. The particle spectra of Lin from which equations (4.1) - (4.4) 
were derived apply to the portion of (4.8) in the range of solid angle 

1 

sin“^120° 
in =- 2TT r de sin 0 
sin" 


and so cannot be used for detefming G(|J.). Supplemental pitch-angle distributions 
indicate a degree of anisotropy that is hard to model analytically, but we 
feel that the distribution, at least over the rise time of the event, is closer 
to that depicted in Figure- 4b. 
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With regard to point (ii), we have already noted above that an integra- 
tion of the flux distributions observed by Lin, assimiing them to be isotropic 
in the forward hemisphere G < tt/ 2, tends to underestimate the density, if in 
fact it is dominated by low-pitch-angle particles. The factor by which the 
density is thus underestimated is difficult to determine precisely, but an 
inspection of Lin’s unpublished date indicates that it may be between one 
and ten. On the other hand, the asstimption noted above that is taken 

to be at its fully evolved value tends to overestimate the density, because 
the peaks we identify in the time-evolving spectra are in fact at slightly 
lower flux levels — for a given energy — than that of the asymptotic spectrum 
(c.f. Figure 3 Lin, Evans, and Fainberg 1973) by a factor which is 

typically about three. Thus, although we do not attempt in practice to correct . 
for either of these errors, their combined effect is to produce an inaccurate 
estimate of the density for any particular event by a factor of order 

unity, which is not important in the context of the type of comparison to 
observations that we shall make below. Moreover, because there is a wide range 
of densities possible in type III exciters in general, the normalizations we 
use are appropriate for illustrating the general physics independent of 
detailed comparisons to particular events. We note also that attempts to 
derive an overall normalization from Lin's data are also subject to uncer- 
tainties by factors of order unity. 

With regard to point (iii), we note first that the schematic depiction 
of f(u) as in Figure 4b includes only the effect of the pitch-angle distribu- 
tion on particles with speed v > v^. The error bars in the flux distributions 
for E^ <E preclude developing a model for dj/dE in this region, although 

the measurements allow the identification of according to equation (4.l). 

Thus, we have merely modelled f(P) in the unstable region by a ramp. 
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The confinement of the ramp to ^ P ^ is consistent with our earlier 
assumptions because the steeply rising portion dominates the exponential 
growth of the linearly unstable spectrum. In addition, it must be remembered 
that observations of the radio emission integrate the emission from a greatly 
extended exciter, which typically has lateral dimension of about 1 A.U. at 
1 A.U. Thus, the effects of fine structure in f(P) are unobservable in 

I 

principle, because such structure would not be spatially uniform. Such fine 
structure, however, has an effect locally on our calculations, because xt 
affects the ten^oral development of the piimp spectrum. In addition, the 
nonlinear processes themselves may affect f(P) in the resonant region. For 
all these reasons, it may easily be appreciated that the model of a ramp for 
f(P) means that we do not treat the interaction of neighboring modes in the 
resonant region correctly in detail, nor can we. 

A final point is that, as can be seen in Sec. Ill, the Landau damping 
rate in the phase-velocity range P contributes to the spectral dynamics. 

In principle, the beam contribution from small-pitch-angle particles with 
P < P^ can modify this rate, weakening Landau damping or even contributing 
slow linear growth to sustain the spectrum at these phase velocities. We 
have assumed that this is a small effect and do not consider it further. 
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V. ENERGY LOSS OF THE EXCITER DURING PROPAGATION 


Because of the disparity hetween the plasma-wave group velocity and 
the velocity of the exciter electrons, local comparisons of the energy 
densities of waves and particles are meaningless. Thus, to evaluate the 
adequacy of the theory from the standpoint of energetics, we must estimate 
the total energy lost in the exciter propagation volume and compare it to 
the total beam energy flux. 

The total energy lost between the injection near the solar surface 
and the point R, say, is given by 


dW(r t) 

AE = J drA(r) j dt — (5.I) 

R t-(r) dt 

o i 

where A(r) is the source area at r, and t^(r), t^(r) are the times at which 
instability at r begins and ends, respectively. For simplicity, equation 
(5«l) implies uniform conditions across the area; i«e, , dW/dt is an average 
over A(r). 

fV 

From our theory we may estimate dW/dt reliably to within an order of 
magnitude on quite sin5)le physical arguments. In the resonant region the beam 
loses the energy necessary to trigger the OXSI, after which this energy is 
quickly scattered into the non-resonant region and the linear instability 
starts again in the resonant region after the ion fluctuations have been 
damped. Because ail the beam energy loss occurs in the resonant region, we 
may consider that it occurs at the steady rate 


dt T 


W V 

T T b 


A ^ 
c 


( 5 . 2 ) 


where, as above, is the OTSI threshold energy density and A , t are 
defined by 


W„ = W e"c = W e^b c 
To o 
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Using the definition P (x^t) = x/ct as in Sec, IV^ we have the transformation 


/ 






dt£(x,t)=^f 

The beam evolution model of Sec. IV gives 




(5.U 


p^(x,tj^(x)) 



(5.5) 


This is not a simple or convenient function of As we shall see below, 

however, most of the beam energy loss occurs at low where T| is largest; 
and there ) tends to values of order 5“ 10* Thus we shall take this 

factor as simply a constant in the calculation below. 

Then equation (5*l) becomes 


AE(r ^ R) 


/ \2 R 

i . 


drA(r)U)^(r) 


mi 


P® (r) W (P )TKP ) 

jQ 10 O 

p 


max/ \ 
P (r) 


P 


(5.6) 


The limits on the integral over p require discussion. The value p”’‘”(r) 

o o 

is the velocity at which the linear instability stops owing to merging of 

the beam distribution into that of the solar wind, so that 5f /&u ^ 0 . The 

linear instability, of course, will be quite weak before the beam evolves 

down to P = so that use of this value for the limit of the integral, 

o o 

together with equation (5*^)j will tend to overestimate AE. The limit 


P 


max 


is the highest value of P for which the beam do’es not evolve through 


the velocity intieirval P, (P )-P in a time less than t (p ), it will turn 

\ D O O CO 


out that the integral in equation (5 ‘6) will not be sensitive to p°*^, but 


is very sensitive to through a factor (p"*^“) where ^ is again the 


rain^ -C 
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spectral index of the beam. Thus, e.g. , for Q = 6, wc will overestimate 


AE 'by 'about two orders of magnitude If we underestimate * P by a factor 

o 

of two. We shall discuss below the estimation of in more detail. 

o 


From Sec. 'IV, we write 


Pv(Pj 


lUP^) = N(P) J p 


b^^o^ ^ P-P. 






dp 


N(P) 


-P 


-G+2 


-P^ CC-i)(C-2) 

0.6 


(5.7) 

i 


where we have approximated p, (P ) £s p and neglected a small term, and 

D O , O 

where 


_ C-1 

N(P) = \(P ^ P)(C-1>(P) 

Furthermore, we shall make the following approximations: 

(i) In equation (5-5)^ we take = const (of order 10, say) and write 


W^ = Wn(r)T^(r) ( 5 . 8 ) 

with W = const; 

(ii) In equation (5*7) we write 


1 


P 



2 

oTE ^ 


(5.9) 


which is good for p^ O.J; 

(iii) Lastly, we take the source area A(r) to be given by S 

A(r) = C2(r)(r-R^)^ , 

If □ is constant, the source is a cone of solid angle H with vertex at the 
solar surface r If 0 is, for ex^ple, piecewise constant over an 


5 ^ 



interval in r, then the source volume is a portion o£ such a cone over this 
interval . 


Defining P = f ®nd using these approximation^ equation (5*6) 


becomes 


(R/R J 


lAE(r < R)N(W,A^,C)J dp(p-l)^n^^''^(p)T^(p)N(P)n(p)[p^(p)] 


■C40.4 


(5.10) 


where 


G(W,A^,C) 


2TTkgR 


(¥)Uj 


W 


a^(C-i)(C-2)(C-oA) 


.22 


W 


where kg is Boltzman's constant^ and the last estimate assumes (v^/ ) 

— 10, AE is expressed in ergs and T (p) in °K. 

To evaluate the integral in equation (5*10), let us assume a model of 
the corona in which there are J partitions such that 

-V 


(5.11) 

2 


and take 


Then 




T(P) = T(P^) 3 


P ^ P P ^ P 

f',2 0 


j = 1,2- . J 


n(p)=fij, ^ Vj 


I ae| g 


fl. n T I p 

jn A* V 


^ -3V./2-T. p ^-C+0.4 

j dp(p-l)^(p/pjj^P ^ ^N(p)[p^ (p)J .J 


(5.12) 



We remark that N(P) is left under the integral because "under this model we 
expect 


%(•=)- Vi) (^) ^ 

whereas below some point P ( i ♦e. ^ P ^ P) we shall have '^(P) 2. Then the 

narraalization of the density changes [cf, equation (4.6)1 such that 


T|(P) == !!{?)(£ 


(0 


-2+V(p) 


( 5 . 15 ) 


^min 


Now^ in order to evaluate equation (5* 12) we need (thus over- 

estimating I Ae| ) by noting that instability will certainly have ceased when 
the total one-dimensional plasma distribution function satisfies 



total 


{It) 


background 



beam 


- 0 


(5-14) 


Because we take (Sf/^P)^ as a ramp. 






(5-15) 


Where the subscript b denotes the beam- Equation (5»l4) should be satisfied 
near the- point where 


*b<^) ® *o<^) (5.16) 

\ 

.For these considerations ^ we may assume that dn the inner corona the supra- 
• thermal' portion of the electron distribution is Maxwellian. 

Thus the ambient distribution is 

-P^/2P ^ 

ft = \ 6 , (5-17) 

where the subscript t -refers, either to the total ambient distribution if 
there* is no distinct suprathermal tail, or to such a, tail if dt exists. 
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Using equations (5*15)-(5* 1?)^ equation becomes 


N(P) ^ 

(pmin)0.6lJ_^prnin)0.4J ^ 

Again- approximating [l - s; 2 , and defining y 

o 

we arrive finally at the expression 



(P 


,mxn 




1 

2 


(^+0.8)6«y 


N 


C-0.4 


N(P) 


= 6^1 C. 


( 5 . 18 ) 


L J 

Equation (5*l8 he solved numerically for in the appropriate parti- 

tion region. 

For the sake of definiteness and illustration^ let us consider a two- 
partition coronal model as described in Table I. This model is in part 
derived after the density model of Riddle (197^)* Using this model, we may 
evaluate equation (5. 12) straight forwardly in terms of at the lower 

bouidary of each region. To determine , we consider that in the inner 
corona the ambient density is Maxwellian with the teii 5 >erature dependence of 
Table I. We use equation (5.15) to scale the beam density back to the inner 
corona,using as parameters at 1 A.U. the t 3 ^ical values 


T 1 = 10"^, c = 5/P = 0.1, ( 5 . 19 ) 

and assuming C! = 1 sr throughout the propagation volume of the beam. It is 
then straightforward to solve equation ( 5 . 18 ) numerically to find 


pmin pmin (fp^g.l) ^ 0.15 

o o 


(5.20) 


Using these values in the integration of equation (5.12) we find the energy 
losses 1 AE 1 in region 1 and 2 , respectively, to be 
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( 5 . 21 ) 


Iae^I 10^\ ergs 

1 Ae^'I ~ 10 ergs 
-k r 

Our calculations indicate W ^ 10 initially uc£. equation (5.29) and 
Sec. II o£ Goldstein, et al . (197^)].' 11^® total energy in the type III 
exciter has been estimated (Lin and Hudson 19T1> Kane and Lin 1972, Lin 

28 -4 

197^) to he of order 10 ergs. Thus, with W — lO” the exciter loses only 
a few percent of its energy. Ninety percent of the energy loss occurs in 
the frequency range below 15 MHz, thus very weak exciters would give rise 
to bursts which would end in this range, as is frequently observed. We 
may conclude in general that our theory is in excellent accord with 
energetic requirements. 

There is an even more striking conclusion to be drawn from this cal- 
culation, however, which affects the general interpretation of observations. 
Recall that the exciter velocity is observationally inferred from the 
frequency drift rate of the peak of the time-intensity profile at fixed 
frequencies. We refer to the velocity thus inferred as the ''nominal 
velocity". It is usually stated that at high frequencies (> I5 MHz) the' 
nominal velocity is ~ c/3 ( cf . Wild and Smerd, 1972), while Evans, Fainberg 
and Stone (1975) found that at low frequencies ^ 1 MHz the nominal velocity 

decreased by a factor of about two between r ~ 10 R and r 1 A.U, 

o 

Note first that the peak intensity at any frequency is reached at or , 
just above the linear beam-plasma instability stops at that frequency; 
this is just the point which we have, calculated above, when 5^ "• P™^”'(r). 
Recalling that our model of the beam evolution gives the result 

0, (P ) 0°'^ 

bo o 


58 



we find B =0.32 for B-. = 0.l6. If we determine P at 1 A.U.^, assuming 

D u O 

that the high-velocity tails of the solar wind electrons are described by 

Maxwellians with T — 1.2 10^ K and a density fraction of — 0.05 (W.C. 

e 

Feldman, private communication), then we find p^ — 0 .OT — 0.2. 

These calculations suggest the following important conclusions: 

(1) The nominal velocity c/3 is not the characteristic velocity to 
which electrons are accelerated in flares or other active phenomena, 
but reflects the dispersive evolution of a particle spectrum which is 
monotone-decreasing when fully evolved. 

(2) The observations of Evans, Fainberg, and Stone (1975) i^ot 
indicate an actual deceleration of the exciter, but rather reflect the 
fact that the solar wind ten^erature decreases slowly with radius, so 
that merging of the beam and ambient distributions occurs at lower 
velocities. 

There is one caveat to be made regarding these conclusions. Although the 

numerical coincidence of the derived 0, with the nominal values is striking, 

□ 

it is of course somewhat dependent on the beam parameters. This dependence is 
weak, however, being manifested only in the term in equation (5.I8). Thus, 

it seems safe for these purposes not only to adopt the typical beam parameters 
(5*19); to neglect consideration of changes in the parameters during 
propagation. To the extent that the model is valid, then, the small charac- 
teristic dispersion of the observed nominal velocity is due to the weak 

0 .6 

dependence of on 6^; . This law was derived in Sec, IV, how- 

ever, only on the basis of observations near 1 A.U., for a few cases of 
beam evolution in which L (the only parameter in applying the beam evolution 
model to a particular event) varied only about 20^. In the calculation of 
this section, however, we have used this law all the'wa:y through the corona. 
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That the results are so benign, to our estimates of energetics, and to the 

correspondence with observations of the nominal velocity, indicates that 

such a weakly dependent -connection between P and p, probably describes 

o D 

the integrated effects of pitch-angle scattering and reabsorption of the 
wave energy by the beam particles during their propagation. This finding 
may be of interest to future theoretical work. 

V. SUMMARY 

We have presented the complete theory of modulational and collapse 
instabilities in three dimensions. As in Paper I, we again find that for 
parameters appropriate to the interplanetary medium, these strong in- 
stabilities are essentially one dimensional. The pump waves excited by the 
linearly unstable electron beam were found to be confined to a narrow pro- 
pagation cone aligned along the ambient magnetic field. The daughter waves 
produced hy the modulational instabilities were also found to be one dimen- 
sional: parallel and anti-parallel to the direction of B. Furthermore, 
the decay instability was shown to have a higher threshold than the OTSI. 
Consequently, the criticisms made by Bardwell and Goldman (19T6) about the 
inadequacy of a one dimensional treatment of modulational instabilities are 
not germane to type III bursts in the interplanetary medium. 

As in our earlier treatments of the initial stages of the OTSI (Papers 
I and II ), spectral transfer of Langmuir waves resonant with the electron 
beam to short wavelength, nonresonant, waves takes place on a time scale that 
is short compared with the time scale for plateau formation and exciter energy 
loss associated with the quasilinear kinetic theory. The present paper has 
extended those results to now include a description of the processes which in 
turn stabilize the OTSI'; enhanced scattering of Langmuir waves by nonthermal 



ion flticttiatioas that have modified the Bohm-Gross dispersion relation of 
the pump and daughter Langmuir waves; and an anomalous resistivity created 
by these same ion fluctuations. 

While the derivation of the equations of spectral transfer which in- 
clude all of the above effects has been given above, the solution of these 
equations using parameters appropriate to solar type III bursts observed i,n 
the interplanetary medium has been deferred to the accompanying paper 
(Goldstein, et al. IS/jQ). in addition to the detailed development of modu- 
ulational instabilities, we have also presented an analytic model that 
describes the energy losses suffered by an electron exciter of tjrpe III 
bursts that is stabilized by the OTSI. Most of the energy loss is found to 
occur above the I5 MHz level. If the electrons have a total energy of some 
10 ergs, then we find that only a few percent of that energy is lost to 
Langmuir waves as the beam propagates out to 1 A.U. Finally, the apparent 
deceleration of the electron exciters that has been observed by Evans, 
et al . (1975) can be interpreted in the context of this model calculation as 
being due to the decrease in solar wind temperature with heliocentric distance. 
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Figure Captions 


Figure 1. 


Figure 2. 


Figure 3* 


Figure 4. 


Schematic plot of a one dimensional electron distribution function 
with a linearly unstable "bump" centered about v^. 

The final state of the distribution function shown in Fig. 1 as 
predicted by the quasi linear kinetic theory of an infinite, 
homogeneous plasma. 

The effect of "oscillation pileup" which occurs when the electron 
beam is injected into a half infinite plasma. Again, the final 
state (b) is that predicted by quasilinear theory. 

Various possibilities for reduced one dimensional distribution 
functions resulting froma simple first order pitch angle aniso- 
tropy (a) proportional to p; or, (b) proportional to p ;n » 1. 
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TABLE I - CORONAL MODEL 


n(P^)(cm’5) f^(P^)(MH^) TjP^)roK) 

Q ^ 

8X10 250 7.5 0.26 2.0x10 

5 X 10^ 15 2.6 0.26 1.6 X 10^ 
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Figure 2 




Figure 3 
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